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Abstract 

Complete eigenvalue spectra of the staggered Dirac operator in quenched 

Ad compact QED are studied on 8 3 x 4 and 8 3 x 6 lattices. We investigate 

the behavior of the nearest-neighbor spacing distribution P{s) as a measure 

of the fluctuation properties of the eigenvalues in the strong coupling and 

the Coulomb phase. In both phases we find agreement with the Wigner 

surmise of the unitary ensemble of random-matrix theory indicating quantum 

chaos. Combining this with previous results on QCD, we conjecture that quite 

generally the non-linear couplings of quantum field theories lead to a chaotic 

behavior of the eigenvalues of the Dirac operator. 
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I. MOTIVATION 



The fluctuation properties of the eigenvalues of the Euclidean lattice QCD Dirac operator 
have attracted much attention in the past few years. In Ref. |i| it was first shown for 
SU(2) lattice gauge theory that certain features of the spectrum of the Dirac operator are 
described by random- matrix theory (RMT). In particular the so-called nearest-neighbor 
spacing distribution P(s), i.e. the distribution of the spacings s of adjacent eigenvalues on 
the "unfolded" scale (see below), agrees with the Wigner surmise of RMT. According to 
the Bohigas-Giannoni-Schmit conjecture ||, quantum systems whose classical counterparts 
are chaotic have a P(s) given by RMT whereas systems whose classical counterparts are 
integrable obey a Poisson distribution P(s) = e~ s . Therefore, the specific form of P(s) 
is often taken as a criterion for "quantum chaos". However, there is no accepted proof of 
the Bohigas-Giannoni-Schmit conjecture yet. The field of quantum chaos is still developing 
and there are many open conceptual problems ||. Applying this conjecture it was recently 
demonstrated that QCD is chaotic, both in the confinement and the quark gluon plasma 
phase 0. 

A number of interesting results have been established for chaotic dynamics in classi- 
cal gauge theories. Lattice gauge theories are chaotic as classical Hamiltonian dynamical 
systems |5]]. Furthermore, it was found that the leading Lyapunov exponent of SU(2) Yang- 
Mills field configurations indicates that configurations corresponding to the deconfinement 
phase are chaotic although they are less chaotic than in the strong coupling phase at finite 
temperature ||. The scaling of the maximal Lyapunov exponent in the classical continuum 
limit was studied in Ref. f7j: It was suggested that Abelian gauge theories behave regularly 
in the continuum limit whereas non-Abelian gauge theories are chaotic in the continuum, 
although the exact scaling relation is still an open problem. Chaos to order transitions were 
observed in a spatially homogeneous SU(2) Yang-Mills-Higgs system and in a spatially ho- 
mogeneous SU(2) Yang-Mills Chern-Simons Higgs system J8]||. In Ref. || a chaos to order 
transition was also seen on the quantum level, i.e. a smooth transition from a Wigner to a 
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Poisson distribution was found. A transition in P(s) from Wigner to Poisson behavior was 
further observed at the metal- insulator transition of the Anderson model [H| . Recently, the 
suppression of the characteristic manifestations of dynamical chaos by quantum fluctuations 
was analyzed in the context of spatially homogeneous scalar electrodynamics |TT| and for a 
+ 1-dimensional space-time iV-component 4 theory in the presence of an external field 



12| . These chaos to order transitions were seen in spatially homogeneous models and not 



for the full classical field theory. The relationship to properties of the quantum field theory 
is an interesting issue. 

Here we focus on the Dirac operator for quenched 4c? compact QED to search for the 
possible existence of a transition from chaotic to regular behavior in Abelian lattice gauge 
theories. In particular, we are interested in the nearest-neighbor spacing distribution of the 
eigenvalues of the Dirac operator across the phase transition from the strong coupling to 
the Coulomb phase. In the strong coupling region Abelian as well as non-Abelian lattice 
gauge theories are in a confined phase ||13|| . For compact QED this means that for couplings 
(3 < f3 c ~ 1.01 the electron is confined. However, when crossing the phase transition the 
conventional Coulomb phase is observed. There are some interesting properties of the two 
phases which can be studied in lattice QED. In the confinement phase the photons form 
massive bound states similar as the gluons bind to glue-balls in lattice QCD. When crossing 
the phase transition a massless photon is found |14[ whereas in lattice QCD the gluon is 
a massive particle in the deconfinement region. U(l) lattice gauge theory contains Dirac 
magnetic monopoles in addition to photons and it was demonstrated via numerical 
simulations that the vacuum in the confined phase is populated by monopole currents which 
become rare in the Coulomb phase [16| . It is an interesting question if the difference between 
the Coulomb phase in QED and the quark-gluon plasma phase in QCD has an influence on 
the level repulsion of the corresponding Dirac spectra. 
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II. ANALYSIS 



We generated gauge field configurations using the standard Wilson plaquette action for 
U(l) gauge theory, 

S g (U 1 )=(3Y,( 1 - cos ®p) > (2- 1 ) 
p 

where Ui = U XfM = exp(iQ XfM ), with Q Xfl G [— ir, ir), are the field variables defined on the links 
/ = (a;,/i). The plaquette angles are Bp = 9 XjAt + Q x +ji,v — ®x+0,v ~ ®x,v We simulated 
8 3 x 4 and 8 3 x 6 lattices at various values of the inverse gauge coupling (3 = 1/e 2 both in the 
strong coupling and the Coulomb phase. Typically we discarded the first 10000 sweeps for 
reaching equilibrium and produced 20 independent configurations separated by 1000 sweeps 
for each (3. Because of the spectral ergodicity property of RMT one can replace ensemble 
averages by spectral averages [17]] if one is only interested in the bulk properties. Thus a 
few independent configurations are sufficient to compute P{s). 

On the lattice the Dirac operator^) = ft + ieJj± for staggered fermions 

1 4 

M XjX > — - rj Xfl (5 x+ fx }X 'U X}IJj — Sx-fax'Ul ij) (2.2) 

is anti-Hermitian so that all eigenvalues are imaginary. For convenience we denote them by 
i\ n and refer to the A n as the eigenvalues in the following. Because of {^,75} = the A n 
occur in pairs of opposite sign. All spectra were checked against the analytical sum rules 

£A„, = and Y, X l = V i ( 2 - 3 ) 

n A„>0 



where V is the lattice volume [IS|. We further checked our spectra by calculating the chiral 
condensate 

(xx) = V- 1 (Y,(^n + m)- 1 } (2.4) 



for m = 0.04 and found agreement with results in the literature |[L9| . 

To construct the nearest-neighbor spacing distribution P(s) from the eigenvalues, one 
has to "unfold" the spectra. This procedure is a local rescaling of the energy scale so that 



the mean level spacing s is equal to unity on the unfolded scale p0| : One first defines the 
staircase function N(E) to be the number of eigenvalues with A < E. This staircase function 
is decomposed into an average part and a fluctuating part, N(E) = N av (E) + Nq(E). The 
smooth average part is extracted by fitting N(E) to a smooth curve, e.g. to a low-order 
Chebyshev polynomial. One then defines the unfolded energies to be x n = A* r av (£' n ). As a 
consequence the sequence {x n } has mean level spacing equal to unity. Ensemble and spectral 
averages are only meaningful after unfolding. Figure [l] shows a typical staircase function for 
(3 = 0.90 (strong coupling phase) and j3 = 1.10 (Coulomb phase) on an 8 3 x 6 lattice. It 
exhibits a decrease of small eigenvalues due to the restoration of chiral symmetry across the 
transition. 



III. RESULTS AND DISCUSSION 

In RMT one has to distinguish between different universality classes which are determined 
by the symmetries of the system. So far the classification for the QED Dirac operator has not 
been done. Our calculations show that in the case of the staggered 4d compact QED Dirac 
matrix the appropriate ensemble is the unitary ensemble. Although from a mathematical 
point of view this is the simplest one, the RMT result for the nearest-neighbor spacing 
distribution is still rather complicated. It can be expressed in terms of so-called prolate 
spheroidal functions, see Ref. where P(s) has also been tabulated. A good approximation 
to P(s) is provided by the Wigner surmise for the unitary ensemble 

P(s) = ^s*e-^. (3.1) 

We have simulated 8 3 x 4 lattices at p = 0, 0.90, 0.95, 1.00, 1.05, 1.10, 1.50 and 8 3 x 6 
lattices at /3 = 0.90, 1.10, 1.50. All results are similar to those selected for the plots. Figure |2| 
shows the nearest-neighbor spacing distribution P(s) for f3 = 0.90 in the confined phase 
averaged over 20 independent configurations on the 8 3 x 6 lattice compared with the Wigner 
surmise for the unitary ensemble of RMT of Eq. ( |3.1| ). Good agreement is found. According 
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to the Bohigas-Giannoni-Schmit conjecture this means the system can be regarded as chaotic 
in the strong coupling region. Figure |3| shows the nearest-neighbor spacing distribution P(s) 
for /3 — 1.10 in the Coulomb phase again averaged over 20 independent configurations and 
compared with the Wigner surmise (|3.1| ). The agreement of the lattice data with the RMT 
predictions is interpreted as a signal that quantum chaos survives the phase transition. We 
find no deviation up to the maximum coupling considered, (3 = 1.50. 

In the strong coupling phase the result holds down to (3 = 0. Therefore, we tend to 
interpret our, as well as previous results in the sense that the disorder of the gauge 

field configurations |5||| is responsible for the chaotic characteristics of the spectrum of 
the Dirac operator. In contrast to that: The free fermion theory is non-chaotic and the 
corresponding nearest-neighbor spacing distribution obeys a Poisson distribution. This is 
illustrated in Fig. [| where P(s) is obtained from the analytical eigenvalues of the free Dirac 
operator on a 53 x 47 x 43 x 41 lattice: 




Here a is the lattice constant, is the number of lattice sites in /x-direction, and n M = 
0, L M — 1. We used an asymmetric lattice with being primes and restricted the range 
to {Lfi — l)/2 instead of L M — 1 in each direction to avoid degeneracies of the free spectrum 

m. 

IV. CONCLUSION 

We have analyzed the nearest-neighbor spacing distribution P(s) of the eigenvalues of 
the Dirac operator in quenched 4d QED on 8 3 x 4 and 8 3 x 6 lattices both in the strong 
coupling region and in the Coulomb phase. In both phases we found excellent agreement 
of the lattice data with the Wigner surmise of the unitary ensemble of RMT. Our results 
evidence that the fermions in U(l) gauge theory show quantum chaos in the confined as 
well as in the Coulomb phase. Dynamical fermions are not expected to affect the Wigner 
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distribution as has been demonstrated for SU(3) ||. In accordance with previous findings 
we conjecture that the eigenvalues of the Dirac operator of interacting quantum field 
theories quite generally reveal chaos due to the disorder of the gauge field configurations. 
The free Dirac operator, in absence of a covariant derivative and minimal gauge coupling, 
exhibits regular behavior. 

It would be interesting to study the relationship between chaos to order transitions in 
classical p|-P, |iT|[T2| and quantum field theories. However, this faces several difficulties: The 
available investigations of classical field theories focus mainly on the gauge sector, whereas 
the numerical methods employed here are only efficient for the fermion sector of quantum 
field theory. A similar accurate determination of the eigenvalue spectrum of the gauge sector 
necessitates to construct the corresponding Fock space and to diagonalize high-dimensional 
matrices which seems to be out of reach for Ad QED/QCD. On the other hand, for the 
classical limit fermion sector studies of chaos have not yet been attempted. 
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FIG. 1. Staircase function N(E) representing the number of positive eigenvalues < E for a 
typical configuration of compact U(l) theory on an 8 3 x 6 lattice in the strong coupling phase 
j3 = 0.90 (solid line) and in the Coulomb phase (5 = 1.10 (dotted line). 
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FIG. 2. Nearest-neighbor spacing distribution P(s) of the Dirac operator for compact U(l) 
theory in the strong coupling phase for (3 = 0.90. The histogram represents the lattice data on an 
8 3 x 6 lattice averaged over 20 independent configurations. The full curve is the Wigner distribution 
of Eq. Q3.1| ) for the unitary ensemble of RMT. 
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FIG. 3. Nearest-neighbor spacing distribution P(s) of the Dirac operator for compact U(l) 
theory in the Coulomb phase for (3 = 1.10. The histogram represents the lattice data on an 8 3 x 6 
lattice averaged over 20 independent configurations. The full curve is the Wigner distribution of 
Eq. ( |3.lD for the unitary ensemble of RMT. For comparison the Poisson distribution P(s) = e~ s is 
also indicated by the dashed line. 
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FIG. 4. Nearest-neighbor spacing distribution P(s) of the analytically calculated eigenvalues 
of Eq. ( |3.2|) for a free Dirac operator on a 53 x 47 x 43 x 41 lattice (histogram) compared with the 
Poisson distribution P(s) = e~ s (solid line). 
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